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I.  Objectives 

Current  research  into  the  development  of  receivers  for  wireless  Code-Division  Multiple-Access 
signals  had  focused  on  systems  where  perfect  knowledge  of  communication  parameters  exists  prior 
to  the  research  funded  by  the  grant  in  question.  It  is  clear  that  such  a  scenario  is  unlikely  as 
most  parameters  will  be  transmitted  over  the  noisy  wireless  channel  or  will  be  estimated.  While 
there  had  been  work  evaluating  the  impact  of  imperfect  knowledge  of  the  necessary  parameters  on 
previously  derived  receivers,  there  was  yet  to  be  research  performed  on  the  development  of  receivers 
that  are  insensitive  to  parameter  estimate  errors. 

The  primary  aim  of  the  research  was  to  investigate  methods  of  introducing  robustness  into 
Direct-Sequence  Code-Division  Multiple- Access  (DS-CDMA)  receivers  that  exploit  both  some  pre¬ 
liminary  knowledge  of  parameter  errors  as  well  as  the  multiple  access  interference  inherent  to  such 
multiuser  systems.  In  addition,  preliminary  research  on  the  novel  area  of  space-time  block  coding 
for  DS-CDMA  systems  was  also  conducted.  Transmit  diversity  (coding  across  multiple  transmit  an¬ 
tennae)  is  another  method  of  increasing  robustness  and  thus  improving  performance  in  DS-CDMA 
wireless  channels.  To  date,  there  appears  to  be  no  work  in  this  area  that  is  specifically  focused 
on  designing  space-time  block  codes  for  DS-CDMA.  There  has  been  work  devoted  to  applying  pre¬ 
viously  proposed  narrowband  codes  to  the  wideband  problem.  By  exploiting  properties  of  the 
DS-CDMA  signals,  we  have  been  able  to  improve  performance. 


II.  Accomplishments 


A.  Robust  MMSE  Receivers 

We  have  developed  an  improved  minimum-mean  squared-error  (MMSE)  based  CDMA  receiver 
which  relies  on  parameter  estimation  techniques  and  a  resulting  probability  density  function  for 
the  estimated  parameters.  We  have  shown  that  the  previously  proposed  MMSE  (assuming  perfect 
parameter  knowledge)  is  most  sensitive  to  imperfections  in  timing  information,  thus  our  focus  has 
been  on  delay-error-insensitive  receivers.  The  delay  error  has  been  modeled  and  using  a  Bayesian 
framework,  the  IMMSE  receiver  developed.  This  receiver  has  been  compared  and  contrasted  to  the 
conventional  MMSE  receiver.  The  probability  of  error  as  well  as  asymptotic  performance  measures 
have  been  computed. 
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Current  research  into  the  development  of  receivers  for  wireless  Code-Division  Multiple-Access  signals  has  focused  on 
systems  where  perfect  knowledge  of  communication  parameters  exists.  It  is  clear  that  such  a  scenario  is  unlikely  as  most 
parameters  will  be  transmitted  over  the  noisy  wireless  channel  or  will  be  estimated.  Past  work  on  a  new  CDMA  receiver 
based  on  tile  minimum-mean  squared-error  (MMSE)  criterion  using  a  Bayesian  framework  has  been  farthered.  In  the  prior 
work  it  was  established  that  the  conventional  MMSE  receiver  is  most  sensitive  to  mismatch  in  time  delay  information.  The 
improved  MMSE  (IMMSE)  receiver  requires  statistical  information  about  the  time  delay  estimators  employed.  In  addition  an 
analysis  of  a  promising  time-delay  estimator  based  on  the  MUSIC  algorithm  was  conducted.  The  new  analysis  is  alternative 
and  further  to  previous  analyses.  Ill  particular,  bias  estimates  are  provided  and  better  approximations  on  the  variance  of  the 
estimators  are  developed.  Recently  the  research  has  focused  on  the  creation  of  timing  error  insensitive  receivers  which 
approximate  maximum-likelihood  sequence  detectors  (MLSD).  Thus  near-optimal  receivers  can  be  achieved  for  mismatched 
asynchronous  channels.  Current  effort  is  directed  towards  a  theoretical  performance  analysis  of  the  modified  MLSDs  and 
towards  the  consideration  of  fading  and  multipath  channels.  Finally,  space-time  block  codes  to  achieve  transmit  diversity 
specifically  for  DS-CDMA  systems  are  obtained  and  analyzed. 
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Figure  1:  AME  of  CMMSE  and  IMMSE  receivers  with  Gaussian  and  uniform  p[y)  in  various 
near-far  situation  (N=31;  6  users;  av  =  v  =  0.2 Tc;  SNR\  =  ^  =  13 dB).  The  IMMSE  receiver 
dominates  in  large  near-far  situation. 


Figure  2:  Cost  ratio  of  direct  inverse  over  the  proposed  mechanism  for  various  ratio  of  code  length 
over  simultaneous  updates. 

The  multi-user  performance  measure  of  interest  in  this  study  was  the  asymptotic  multi-user  efficiency 
(AME)  The  AME  describes  the  performance  degradation  of  a  receiver  due  to  the  multiple-access 
interference  as  the  additive  channel  noise  diminishes  for  a  fixed  set  of  received  powers  for  all  users. 
Based  on  the  AME  definition  in,  Figure  1  shows  the  performance  of  the  conventional  MMSE 
(CMMSE)  and  IMMSE  receivers  with  both  Gaussian  and  uniform  distributions  for  v.  The  advan¬ 
tage  of  the  IMMSE  receivers  for  large  near-far  ratios  is  well  demonstrated. 

A  method  for  dynamically  updating  the  IMMSE  receiver  as  channel  statistics  vary  has  been 
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proposed.  Figure  2  shows  the  computational  cost  ratio  as  a  function  of  the  code  length  N  over 
the  number  of 'simultaneous  updates.  For  example,  when  a  system  with  N=255  has  a  new  user 
coming  in,  the  proposed  update  scheme  requires  40  times  less  computation  than  direct  inversion. 
The  proposed  mechanism  has  also  been  extended  to  accommodate  a  delay  estimator  with  varying 
statistics. 

B.  Time  Delay  Estimator  Analysis 

Motivated  by  prior  work  which  developed  an  improved  minimum-mean  squared-error  (MMSE) 
based  CDMA  receiver  which  is  less  sensitive  to  errors  in  timing  information  than  the  conventional 
MMSE  receiver,  the  analysis  of  timing  estimators  for  CDMA  channels  has  been  considered.  In 
particular,  the  probability  density  functions  of  the  modified  MUSIC  delay  estimator  has  been 
performed.  An  analysis  that  is  alternative  and  further  to  that  previously  performed  has  been 
conducted.  In  particular,  the  bias  of  the  modified  MUSIC  estimator  for  the  CDMA  channel  has 
been  evaluated  (this  quantity  had  not  been  previously  studied)  and  improved  estimates  of  the 
variance  of  the  modified  MUSIC  estimator  were  also  developed. 

In  general,  the  exact  properties  of  the  time  delay  estimates  with  finite  sample  size  or  low  signal- 
to-noise  ratio  (SNR)  are  intractable.  In  [6],  Parkvall  provided  an  asymptotic  performance  analysis 
of  the  modified  MUSIC  estimator  [5].  In  the  current  work,  an  alternative  perturbation  analysis' of 
the  second  order  statistics  was  pursued  which  also  led  to  a  first  order  approximation.  These  new 
results  are  provided  in  our  work.  It  is  noted  that  the  analysis  in  [6]  did  not  consider  the  estimation 
of  the  bias. 

Figure  3  shows  the  analytical  and  simulated  results  of  the  MUSIC  estimator  with  10  users. 
The  interfering  users  are  10  dB  stronger  than  the  interested  user.  The  signal-to-noise  ratio  (SNR), 
defined  as  symbol  energy  over  the  two-sided  power  spectral  density  of  noise  No/2,  is  15  dB.  The  sim¬ 
ulation  results  are  drawn  from  4000  Monte-Carlo  runs.  The  symbols  of  each  user  are  spread  by  Gold 
codes  of  length  N=Z1.  The  delay  set  is  (25.40, 22.90, 26.95, 9.30, 8.28, 16.21, 19.67, 22.48, 29.05, 11.31) 
chips.  It  is  clear  that  the  new  approximations  are  quite  close  to  the  simulated  results  and  further¬ 
more  are  closer  to  the  simulation  data  than  the  approximations  of  [6].  Figure  3  (c)-(d)  shows  the 
same  scheme  as  in  Figure  3  (a)-(b)  with  the  energy  of  each  interfering  user  20  dB  stronger  than  the 
interested  user.  The  new  approximations  are  perceivably  better  than  that  of  [6]  in  all  situations. 
It  is  observed  that  both  the  estimator  performance  and  its  approximation  are  fairly  insensitive  to 
the  energy  variation.  However,  they  are  sensitive  to  the  change  of  code  sequences  and  delays.  The 
operation  of  the  IMMSE  receiver  has  been  studied  using  a  modified  MUSIC  delay  estimator  as  a 
front  end.  The  results  have  been  very  promising  and  further  analysis  is  being  undertaken  in  order 
to  quantify  the  performance  gains  of  this  integrated  system  over  other  systems. 

C.  Robust  Sequence  Detectors 

The  initial  study  on  MMSE  receivers  enabled  the  consideration  of  optimal  multi-user  detector 
for  asynchronous  channels.  That  is,  a  maximum-likelihood  sequence  detector  (MLSD)  have  been 
modified  to  reduce  sensitivity  to  mismatch  in  timing  errors. 

A  modified  MLSD  has  been  designed  and  analyzed.  Furthermore,  a  recursive  implementation 
allowing  for  use  in  practical  systems  has  also  been  derived.  Performance  of  the  original  MSLD  and 
robust  MLSD  algorithms  has  been  evaluated  via  simulation  and  the  determination  of  upper  and 
lower  bounds  on  the  bit  error  rate.  This  analysis  has  occurred  for  both  additive  white  Gaussian 
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Figure  3:  Approximation  and  simulation  of  estimator  performance  with  10  users:  The  interfering 
users  are  10  dB  stronger  than  the  interested  user  in  Figures  (a)-(b)  and  20dB  in  Figures  (c)-(d). 
Dashed  line  with  o  is  simulation;  solid  line  is  the  new  approximation;  dashed-dotted  line  and  dotted 
line  are  the  approximation  of  mean-squared  error  and  asymptotic  Cramer-Rao  bound  by  Parkvallf], 
respectively. 

noise  channels  as  well  as  flat  fading  channels.  In  the  flat  fading  channel  we  see  a  degradation  in 
performance;  this  was  to  be  expected  as  there  is  no  form  of  diversity  in  the  system.  However, 
we  anticipate  that  in  multipath  (frequency-selective)  channels  or  in  systems  with  multiple  sensors, 
performance  would  improve  significantly  for  both  the  original  MLSD  algorithm  as  well  as  the 
robustified  MLSD  algorithm. 

Efforts  towards  reducing  the  impact  have  been  previously  considered  in  [1,  2,  3,  4].  Most, 
if  not  all,  of  the  previous  works  considered  single  or  multi-shot  type  receivers  where  edge  effects, 
defined  as  a  partial  observation  of  the  data  at  the  edge  of  the  observation  window,  contribute  to  the 
degradation  of  detection  performance.  In  pursuit  of  near-optimum  performance  with  complexity 
linear  in  the  sequence  length,  we  propose  a  modified  maximum-likelihood  sequence  detector  (MLSD) 
which  combats  the  effects  of  imperfect  synchronization. 

The  performance  of  the  O-MLSD,  M-MLSD  and  RM-MLSD  is  presented  by  simulated  BER 
and  analytical  bounds  versus  near-far  ratio  (NFR),  SNR,  and  fixed  realization  of  the  delay  errors 
for  two  active  users.  The  m-sequence  of  length  15  with  two  different  phases  is  applied  for  the  chip 
sequence.  The  number  of  sample  for  each  user  in  each  chip  interval  is  3,  i.e.,  qc h  =  3.  The  decision 
depth  of  the  simulated  MLSDs  is  five  times  the  channel  memory.  In  each  figure,  the  solid  lines  are 
simulation  results  for  the  O-MLSD  and  the  RM-MLSD.  Their  corresponding  bounds  are  presented 
by  the  dash-dotted  lines  and  the  dotted  lines,  respectively.  To  quantify  the  performance  loss  by 
the  recursive  implementation  of  the  modified  ML  criterion,  newly  derived  bounds  for  the  optimal 
ML  criterion  are  presented  by  dashed  lines. 

Three  receivers  were  studied:  the  original  MLSD  algorithm  which  ignores  parameter  mismatch 
(O-MLSD);  a  modified  MLSD  algorithm  (M-MLSD)  which  has  high  complexity,  but  incorporates 
the  effects  of  the  mismatched  channel;  and  the  RM-MSLD  (recursive  M-MLSD)  which  incurs  further 
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Figure  4:  BER  and  its  bounds  vs.  SNR  for  the  O-MLSD  and  the  (R)M-MLSD.  The  parameters  are: 
2  active  users;  N=15;  NFR=20  dB;  aAr  =  0.6TC;  |Ar|  =  0.2TC;  decision  depth=10;  and  qch  =  3. 

approximation  to  provide  a  recursive  implementation.  In  Figure  4,  the  BER  is  demonstrated  as 
a  function  of  SNR  with  NFR=20  dB  and  fixed  delay  errors  of  0.2TC.  The  loss  due  to  recursive 
implementation  is  limited,  especially  in  the  low  SNR  regime  where  the  receiver  is  mostly  challenged. 

D.  Space-Time  Block  Coding  for  DS-CDMA 

Transmit  diversity  methods  have  proven  effective  for  combating  fading  in  wireless  communica¬ 
tion  systems  [7,  9,  11,  10].  In  recent  work,  the  performance  criteria  and  design  of  space-time  block 
codes  have  been  addressed  for  single-user  systems.  There  has  been  some  preliminary  work  to  gen¬ 
eralize  the  idea  of  space-time  coding  to  DS-CDMA  systems  [9,  8].  These  preliminary  investigations 
have  typically  applied  narrowband  space-time  coding  ideas  to  CDMA  systems  and  do  explicitly 
take  advantage  of  the  properties  of  DS-CDMA. 

We  have  developed  a  general  model  for  space-time  modulation  in  CDMA  systems.  From  the 
Chernoff  bound  of  the  probability  of  decoding  error,  the  performance  criteria  are  derived  and  the 
design  of  optimal  space-time  block  codes  (STBC)  is  discussed.  The  presence  of  spreading  codes  in 
the  CDMA  problem  yields  interesting  differences  in  code  designs  and  code  metrics  relative  to  the 
single-user  narrowband  case  [7,  10].  Hence  the  approach  for  CDMA  code  design  is  quite  different 
from  the  narrowband  case. 

Optimal  codes  and  equivalent  classes  are  determined.  Guidelines  for  optimal  code  design  are 
discussed.-  We  have  shown  that  the  code  design  is  completely  decoupled  from  user  to  user,  and 
that  full  diversity  gain  is  easily  achieved  and  coding  gain  is  a  monotonically  decreasing  function 
of  spreading  code  correlation  |p|.  Optimal  codes  are  presented  for  some  simple  cases,  it  is  found 
that  non-unitary  codes  usually  outperform  unitary  codes  in  DS-CDMA  systems.  A  systematic  way 
to  generate  optimal  codes  is  still  an  open  question.  Performance  analysis  by  means  of  numerical 
integration  has  also  been  accomplished.  By  examining  the  union  bound  on  the  probability  of  error 
of  the  ML  decoder  (see  Figure  5),  and  by  simulation,  we  can  see  that,  Class  1  will  be  better  than 
Class  2  when  |p|  is  less  than  about  0.6.  For  CDMA  systems,  |p|  is  usually  designed  to  be  small,  so 
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Figure  5:  Union  bound  vs  p. 


we  expect  that,  in  general,  Class  1  is  better  for  CDMA  systems  than  Class  2. 

Figure  6,  the  numerical  integration  results  are  compared  with  simulation  results  for  BPSK  Class 
1  codes.  The  results  match  very  well. 
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Abstract  —  Space-time  block  code  design  and 
decoder  design  are  addressed  for  Code-Division 
Multiple-Access  (CDMA)  systems.  Optimal  code  de¬ 
signs  are  found  by  optimizing  the  ChernofF  bound  of 
the  probability  of  decoding  error.  From  this,  the  di¬ 
versity  gain  and  the  coding  gain  are  determined  for 
the  CDMA  scenario.  The  resultant  optimal  code  de¬ 
signs  sire  classified  and  analyzed.  Both  optimal  and 
moderate  complexity  suboptimal  decoding  algorithms 
are  proposed  and  evaluated. 

I.  Introduction  &  System  Model 

Transmit  diversity  methods  have  proven  effective  for  com¬ 
batting  fading  in  wireless  communication  systems.  In  this  pa¬ 
per,  we  focus  on  determining  space-time  block  coding  meth¬ 
ods  [1]  for  Code-Division  Multiple-Access  systems.  Due  to 
the  assumption  of  independent  fading  for  different  users ;  the 
multiuser  coding  problem  decouples  to  multiple  single-user  cod¬ 
ing  problems.  The  presence  of  spreading  codes  in  the  CDMA 
problem  yields  interesting  differences  in  code  designs  and  code 
metrics  relative  to  the  single-user  narrowband  case  [4]. 

The  up-link  between  K  users  and  one  base  station  is  con¬ 
sidered.  For  each  user,  kc  information  bits  axe  mapped  to  one 
of  2kc  Space-Time  Block  Codes  (STBC),  D .  The  codeword 
D  is  a  matrix  of  dimension  It  x  M;  where  M  is  the  num¬ 
ber  of  transmit  antennae  and  Lt  is  the  duration,  in  symbol 
intervals,  of  the  code.  It  is  assumed  that  the  base  station 
has  N  receive  antennae.  We  further  assume  that:  the  fading 
processes  associated  with  each  transmit  antenna  are  indepen¬ 
dent;  the  channel  is  constant  over  the  duration  of  the  block 
code  (quasi-static)  and  is  known  perfectly;  and  that  the  trans¬ 
mission  is  synchronous.  For  the  system  under  consideration, 
a  different  spreading  code  is  employed  for  each  transmit  an¬ 
tenna  and  it  is  assumed  that  the  receiver  has  full  knowledge 
of  these  spreading  codes. 

II.  Performance  Criteria  &  Code  Designs 

It  can  be  shown  that  optimizing  the  upper  bound  on  prob¬ 
ability  of  decoding  error  yields  two  criteria  for  space-time 
block  code  design.  Performance  is  determined  by  a  key  ma¬ 
trix  4>  =  A DH  AD  ©  R  lt  where  A D  is  a  codeword  difference 
matrix  and  R  is  the  spreading  code  cross-correlation  matrix. 
The  resultant  design  “metrics”  are: 

diversity  gain  A h  —  where  rmin  is  the  minimum 

rank  of 

coding  gain  Ap  =  (IInmin  is  the  smallest  product  of 
all  the  non-zero  eigenvalues  of  $. 

These  “metrics”  are  analogous  to  those  obtained  in  [1]  for 
the  narrowband  case,  but  due  to  the  presence  of  the  cross- 
correlation  matrix  R ,  some  new  features  appear  in  the  re¬ 
sulting  optimal  codes.  The  goal  of  code  design  is  to  find  2kc 

1Here  O  represents  Schur  product. 


distinct  STBCs  such  that  rm»n  is  maximized  and  given  this 
rank,  that  Ap  is  also  maximized.  Codes  satisfying  these  con¬ 
ditions  are  deemed  optimal  The  following  two  propositions 
can  be  proved  regarding  diversity  gain  and  coding  gain: 

PROPOSITION  1  If  AD  has  no  zero  columns  and  if  R  is 
positive  definite,  full  diversity  gain  is  always  achieved. 

PROPOSITION  2  If  Rij  =  p  for  i  ^  j,  the  coding  gain  is 
a  monotonically  decreasing  function  of  p . 

Thus  for  the  CDMA  case,  we  focus  on  maximizing  Ap.  An 
interesting  observation  is  that  non-unitary  code  usually  out- 
performes  unitary  code.  Consider  the  optimal  code  sets  for 
BPSK  modulation.  We  discuss  the  case  of  kc  =  2,  M  =  2 
and  Lt  =  2.  The  resulting  optimal  codes  can  be  partitioned 
into  three  equivalence  classes.  Each  element  of  the  equiva¬ 
lence  class  can  be  transformed  into  another  element  via  sim¬ 
ple  isometries.  Each  class  is  geometrically  uniform.  Class 
1  and  2  are  optimal  for  all  \p\  <  1  while  Class  3  is  opti¬ 
mal  for  \p\  <  V3/2.  Class  2  is  essentially  Alamouti’s  or¬ 
thogonal  code  set  [3]  while  Class  1  is  non-unitary  [4].  In¬ 
terestingly,  for  QPSK  modulation,  we  can  find  an  optimal 
non-unitaxy  code  set  which  outperforms  all  unitary  codes 
for  all  \p\  <  1.  The  optimal  codes  axe  tabulated  below. 


symbol 

Class 

D1 

D2 

D3 

D4 

BPSK 

1 

1  1 

1  1 

(-•  b-> 

1  1 
W  t~* 

-1  1 

1  -1 

-1  -1 

-1  1 

BPSK 

2 

1  1 

1  -1 

1  -1 

-1  -1 

-1  1 

1  1 

-1  -1 

-1  1 

BPSK 

3 

1  1 

1  1 

1  -1 

-1  1 

-1  1 

1  -1 

-1  -1 

-1  -1 

QPSK' 

1  1 

1  i 

1  -1 
-1  i 

-1  i 
-i  1 

-i  -i 

i  -i 

III.  Decoding  Algorithms 
Three  types  of  decoders  are  considered:  the  optimal 
maxmimum-likelihood  (ML)  decoder,  a  joint  multiuser  mini¬ 
mum  mean-squared  error  decoder  and  a  combined  interference 
cancellation/ ML  decoder.  These  algorithms  perform  as  pre¬ 
dicted  with  the  ML  decoder  offering  the  best  performance  at 
the  expense  of  computational  complexity.  The  two  suboptimal 
algorithms  offer  solid  performance  with  reduced  complexity. 
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Abstract —  Most  receiver  designs  for  asynchronous  direct- 
sequence  code-division  multiple-access  (DS-CDMA)  systems 
exploit  timing  information  to  simultaneously  detect  the  desired 
signal  and  suppress  the  interference.  In  this  paper,  a  previously 
proposed  MUSIC  delay  estimation  algorithm  is  considered  which 
requires  no  initial  information  and  no  exhaustive  search  of  the 
delays.  Based  on  a  Taylor  series  expansion,  approximations  of 
the  first  and  second  moments  of  the  delay  estimation  error  are 
derived  for  this  MUSIC  algorithm.  The  analysis  is  alternative 
and  further  to  that  previously  performed. 

Index  Terms —  Channel  estimation,  code-division  multiple- 
access,  multiaccess  communication. 


I.  Introduction 

A  SIGNIFICANT  number  of  the  previously  proposed  re¬ 
ceiver  designs  for  multiuser  receivers  in  direct-sequence 
code-division  multiple-access  (DS-CDMA)  environments  as¬ 
sumed  that  any  necessary  communication  parameters  were 
perfectly  known  (see,  e.g.,  [1]  and  [2]).  This  parameter  knowl¬ 
edge  was  exploited  to  successfully  suppress  multiple  access 
interference  (MAI)  as  well  as  detect  the  signal  of  the  user  of 
interest.  Of  the  various  relevant  parameter  estimation  problems 
in  a  multiuser  context,  acquisition  of  the  signals  is  one  of 
the  most  challenging.  Thus,  there  has  been  considerable  effort 
devoted  to  the  development  of  delay  estimators  for  DS- 
CDMA  systems.  The  conventional  sliding  correlator  suffers  in 
a  near-far  scenario  when  the  signal  power  of  the  interfering 
users  overwhelms  that  of  the  desired  user  [3].  Other  parameter 
estimation  algorithms  either  assume  initial  information  [4]  or 
require  long  observation  periods  [5].  Our  analysis  is  focused 
on  the  modified  MUSIC  delay  estimator  proposed  in  [3],  which 
requires  no  side  information  about  the  delays.  This  estimator 
is  of  particular  interest  due  to  its  good  performance  with  a 
moderate  number  of  observations.  In  this  paper,  the  study  of 
the  modified  MUSIC  delay  estimator  is  motivated  by  the  need 
for  estimator  statistics  to  drive  the  construction  of  the  delay- 
error  insensitive  minimum-mean  squared-error  (MMSE)  based 
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receivers  proposed  in  [6].  To  properly  construct  the  receivers 
in  [6],  both  bias  and  variance  information  about  the  delay 
estimators  are  necessary. 

In  general,  the  exact  properties  of  the  estimates  with  finite 
sample  size  or  low  signal-to-noise  ratio  (SNR)  are  intractable. 
In  [7],  Parkvall  provided  an  asymptotic  performance  analysis 
of  the  modified  MUSIC  estimator.  In  this  paper,  we  pursue  an 
alternative  perturbation  analysis  of  the  second-order  statistics, 
which  also  leads  to  a  first-order  approximation.  The  analysis 
in  [7]  does  not  provide  an  approximation  for  the  bias  of  the 
estimator,  whereas  the  approach  considered  herein  lends  itself 
to  the  approximation  of  the  bias. 

This  paper  is  organized  as  follows.  In  Section  II,  the 
system  model  and  assumptions  are  provided.  In  addition,  the 
modified  MUSIC  algorithm  is  reviewed.  In  Section  III,  the 
new  performance  analysis  of  the  modified  MUSIC  algorithm 
is  provided.  The  bias  and  mean-squared  error  estimates,  are 
compared  to  previous  estimates  and  simulation  results  in 
Section  IV.  Final  conclusions  are  found  in  Section  V. 

II.  System  Model  and  Modified  MUSIC  Algorithm 

The  system  model  and  the  MUSIC  algorithm  summarized 
in  this  section  are  adopted  from  [3].  In  a  /f-user  asynchronous 
DS-CDMA  system,  each  user  employs  a  spreading  waveform 
bk(t)  with  unit  energy  to  form  the  baseband  signal.  These 
spreading  waveforms  are  formed  by  modulating  rectangular 
chip  waveforms  with  binary  pseudonoise  sequences  [8].  The 
total  received  signal  is  the  sum  of  the  signals  from  all  active 
users  plus  white  Gaussian  noise.  The  baseband  signal  can  be 
represented  as 

K  oo 

r(t)  =  vT^cos ekdk{m)bk{t  -  mT  -  rk)  +  n(t) 

k=l  m——oo 

(D 

where  y/Pk  is  the  amplitude.  Ok  is  the  phase,  r*  is  the  time 
delay,  and  the  data  dk(m)  €  {+1,  —1},  for  BPSK  modulation. 
Without  loss  of  generality,  it  is  assumed  r*.  >  0.  We  assume 
that  the  data  <4(m)  are  identically  distributed  and  independent 
for  all  k  (users)  and  m  (time).  Furthermore,  the  data  and  the 
noise  processes  are  also  independent.  The  received  signal  is 
matched  filtered  with  the  chip  pulse  shape  without  any  timing 
alignment  with  the  delay  Tk-  The  matched  filter  output  is 
sampled  at  the  chip-rate  Tc.  The  samples  corresponding  to 
the  mth  frame  of  the  received  signal  from  the  k\h  user  are 
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given  by 


where  the  received  amplitude  is  fik  =  x/R-cos  Ok-  The  data 
frames  are  of  length  N,  which  correspond  to  the  number  of 
chips  in  the  spreading  waveform.  The  vectors  a2k-i  and  a2k 
represent  the  shifted  signal  sequences, 

a2t_,  =  +  (i  -  <=* 

+  (i  -  <4> 

where  pk  and  6k  denote  the  integer  and  fractional  part  of 
the  delay,  respectively,  i.e.,  rk  =  pkTc  +  Skj  where  pk  is  a 
nonnegative  integer  and  6k  £  [0,  Tc).  The  length  N  vector 
c*  is  the  spreading  vector  of  user  k .  The  permutation  matrix 
Dip)  is  defined  in  block  form  as 

D«=[?  In~A  (5) 

[sip  0 

where  Iat  is  the  identity  matrix  of  dimension  N.  The  vectors 
a2fc_i  and  a2fc  are  modulated  by  different  combinations  of  two 
adjacent  binary  symbols,  z2fc-i  and  z2fc>  respectively, 

z2k_i(m)  =  -(dfe(m)  +  dk{m  -  1)) 

Z2k(m)  =  ^(dk(m)  -  dk{m  -  1)).  (6) 

The  autocorrelation  matrix  of  the  received  signal  r(m), 
where  r(m)  =  Sfcr*(m)  +  n,  is  given  by  R  = 
E{r(m)rH(m)}.  The  superscript  “J?”  denotes  the  conjugate 
transpose.  For  N  >  2  K,  the  matrix  R  can  be  decomposed 
into  nontrivial  signal  and  noise  subspaces.  We  assume  that  the 
signal  vectors  a<,  t  =  1,  •••,  2 K,  are  linearly  independent, 
so  the  signal  subspace  has  dimension  2 K.  In  practice,  the 
autocorrelation  matrix  must  be  estimated.  Given  M  frames 
of  observation  data,  an  unbiased  estimate  of  the  correlation 
matrix  and  its  eigenvalue  decomposition  (EVD)  are  given  by 

R  =  Y"  r(m)rH(m)  =  E5A,Ef  +  E„A„E^  (7) 

m—\ 

where  the  columns  of  E,  and  En  are  the  estimated  eigenvec¬ 
tors  of  the  signal  and  noise  space,  respectively;  A,  and  A„ 
are  diagonal  matrices  of  the  corresponding  eigenvalues  on  the 
diagonal.  In  a  K -user  system  with  N  >  2 K,  E,  is  N  x  2 K 
and  E„  is  N  x  (N  —  2 K).  Their  column  vectors  are  denoted 
by  E,  =  [e,,  •  •  • ,  e2jf]  and  E„  =  [e2/r+i»  •  •  • ,  eAr]- 
In  the  following,  the  user  of  interest  is  labeled  as  user  1 
and  the  subscript  Jfc  =  1  is  dropped.  According  to  the  MUSIC 
algorithm  (e.g.,  see  [9]),  the  delay  estimate  is  determined  such 
that  the  projection  of  the  signal  vectors  onto  the  noise  space  is 
minimized.  However,  the  minimization  of  the  corresponding 
cost  function  J(r)  is  problematic  due  to  the  nature  of  the 
rectangular  chip  waveforms  employed.  To  be  specific,  the  first 
derivative  of  J(r)  with  respect  to  r,  denoted  by  J (t ),  is 
not  continuous  at  the  chip  boundary.  A  practical  approach 


is  to  separate  the  estimation  problem  into  N  subestimation 
problems  [3].  In  each  of  the  problems,  the  delay  estimate 
is  constrained  to  be  within  one  chip  interval,  where  J(t )  is 
analytic,  so  the  MUSIC  algorithm  may  be  applied  directly.  As 
asynchronous  transmission  is  considered,  the  received  signal 
from  the  first  user  is  a  linear  combination  of  two  vectors, 
ai  and  a2.  Hence,  the  generalized  cost  function  for  the  pth 
subestimator  is  a  weighted  squared  norm  of  the  2-D  manifold 
defined  by  Span(ai,  a2)  onto  the  noise  space  E„  [3] 

Jp(Tp)  =  tr|up(rp)A^ (rp)EnEn  Ap(t7,)  j> 

=  tr{Up(rp)A"(rp)(I  -  E.E? )Ap(rp)}  (8) 

where  the  weighting  matrix  Up  is  a  2  x  2  real  matrix,  and 
A p  =  [ai a2]  is  Nx2  with  column  vectors  ai  and  a2  defined 
in  (3)  and  (4).  The  set  Ip  defines  the  interval  [(p  - 1  )TC,  pTc ], 
and  tr{  }  denotes  the  trace  operator.  The  resulting  overall 
estimate  is  given  by 

f  =  arg  min  Jp( fp) 

fi  ■■■  f/^ 

where 

fp  =  arg  min  Jp(rp).  (9) 

Tp€/p 

That  is,  the  estimate  is  selected  from  one  of  the  N  subesti- 
mates.  For  r  =  p0Tc  +  6,  fPo  is  a  consistent  estimate  of  r 
[9],  thus,  f  is  also  a  consistent  estimate  of  r  provided  that 
the  minimum  of  J  is  unique.  However,  for  a  finite  number  of 
samples,  fp  is  not  assured  to  be  the  consistent  estimate  fPo. 
We  shall  refer  to  the  inconsistent  estimates  fp  (p  #  po)  as 
“outliers.”  That  is,  these  outliers  correspond  to  the  incorrect 
chip  interval  altogether.  Since  the  statistical  properties  of  the 
outliers  are  intractable,  the  analytical  performance  provided  in 
Section  HI  excludes  the  outliers  as  done  in  [3]  and  [7]. 

m.  Performance  Analysis  of  the  MUSIC  Estimator 

In  the  remainder  of  this  paper,  we  only  consider  the  proper¬ 
ties  of  the  consistent  estimate  fp,  where  p  =  po,  and  drop 
the  subscript  p.  In  this  section,  we  present  a  performance 
analysis  of  the  modified  MUSIC  estimator.  Our  objective  is 
to  develop  estimates  for  the  first  and  second  moments  of 
the  estimation  error.  This  analysis  is  alternative  and  further 
to  that  considered  in  [7].  The  approach  is  similar  to  [10], 
which  follows  a  perturbation  technique  based  on  a  Taylor 
series  expansion. 

The  Taylor  series  expansion  of  j  about  the  point  (r,  Ea) 
is  given  by 

2  K 

j(f,  Es)  «  j(r,  Ea)  +  J(T,  Ea)A r  +  2Re  £  gf  Ae< 

(10) 

where  At  =  t  —  t,  and  Aej  =  e*  —  ei,  i  =  1,  •  •  • ,  2 K.  The 
vector  g i  is  the  gradient  of  j(r,  Es)  along  e; 

g i  =  VeJ(r,  Ea) 

=  -  {a*UAT  +  A*UAr  +  a*UAT)e*.  (11) 

The  details  of  the  derivation  and  appropriate  definitions  are 
provided  in  the  Appendix.  The  superscripts  *  and  T  denote  the 
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conjugate  and  the  transpose,  respectively.  Given  the  estimated 
signal  space  Es,  the  estimate  of  r  by  the  MUSIC  algorithm 
satisfies 

f  =  arg.  min  J(r,  Es).  (12) 

TG[(p-i)reipre] 

In  the  following,  we  assume  that  due  to  the  definition  of  f 
in  (12),  it  is  necessarily  a  stationary  point  of  J(r,  Es).  Thus 
j(f,  E5)  =  0.  Furthermore,  j(r,  Es)  =  0  where  r  is  the 
true  delay  value  and  E8  is  the  true  signal  subspace.  Thus,  by 
setting  (10)  to  zero,  we  obtain 

2  K 

-2Re£gfAei 

Arw - - .  (13) 

J(r ,  Es) 

The  second  moment  of  At  can  then  be  obtained  by  standard 
manipulation  of  (13) 

2  K  2  K 

•  {gf£’[AeiAeJ]gi  +  gf£[AeiAef]g;}  (14) 

where  E[Aet  AeJ]  and  E[AeiAef  ]  can  be  approximated  by 
two  functions  of  the  eigenvalues  (A j)  and ‘eigenvectors  (e:)  of 
the  true  autocorrelation  matrix  R  [11] 


£[AeiAef 


A  _  V  Ak 
d'j  M  "  (Aj  -  Afc) 

kjii 

.ri  k..\  ~^i 


2e*ek> 


£[Ae,AeJ]  a  (1  -  %)jj^T^2eiei  ' ■  <16> 

The  Kronecker  delta  function  is  denoted  by  Sij. 

Bias  of  the  modified  MUSIC  delay  estimator  of  [3]  has 
not  been  previously  addressed  in  the  literature;  though  bias 
analysis  has  been  considered  for  the  original  MUSIC  estimator 
for  direction  of  arrival  [10],  as  it  is  an  important  performance 
characteristic  for  limited  sample  size  or  low  SNR.  In  addition, 
such  information  is  necessary  for  the  construction  of  the  delay 
error  insensitive  estimators  of  [6].  To  fully  utilize  the  statistical 
properties  of  the  eigenvectors  up  to  the  second  order,  we 
expand  j  by  a  second-order  Taylor  series  expansion  about 
the  point  (r,  es) 

2  K 

j( f,  Eg)  w  J(t,  Eg)  +  J(t,  E,)Ar  +  2Re  ^gfA^ 


+  \  j  (r,  E„)  At2  +  2Rc  ^  [qf  Ae;Ar] 


+  Re tr  ^  HijAejAef  +  H{, jAejAefj 

(17) 

where  the  vector  q*,  calculated  in  the  Appendix,  is  the  gradient 
of  J  along  e* 

q»  =  Vei  J(r,  E5) 

=  -A*UAt  +  2(a*UAt  +  A*UAt  +  A*UaT)  e*. 


The  Hessian  matrices  are  derived  in  the  Appendix  as  well, 

H ij(k,  l) 

a2  j(r,  Eg) 

de*(k)dej(l) 

=  f  -  [(A*UAr  +  A‘UAt  +  A*UAt)]  {  ,  i  =  j 

1 0,  i.^j 

(19) 

H ij(k,  l) 

d2i(T,Eg) 

dei(k)dej{l ) 

=  0.  (20) 

By  setting  (17)  equal  to  zero,  we  obtain  the  approximation 
of  i?[Ar] 

2  k 

E[ At)  «  ..  — [—  Y)  Re{2gf  E[Aei]  +  2qf  E{Aei At] 
J(t,  Eg)  ££ 

+  tr{HME[AeiAef)}}  +  \  J(r,  Es)E[At2} 

(21) 


where  the  bias  of  the  eigenvectors  can  be  approximated  by  [1 1] 

<22) 


The  approximation  of  £[AeiAr]  is  given  by  premultiplying 
(13)  with  Ae;  before  taking  the  expectation  [10] 


E[AeiAr] 


-1 

J(r,  Eg) 


E 


■  {£[aeiAef  ]g;  +  £[Ae,Ae[]g,}.  (23) 


Finally,  the  bias  is  obtained  by  substituting  (14),  (15),  (22), 
and  (23)  into  (21). 


IV.  Numerical  Results 

Figs.  1  and  2  show  the  analytical  and  simulated  results 
of  the  MUSIC  estimator  with  five  and  ten  users.  The  SNR, 
defined  as  the  symbol  energy  over  the  two-sided  power  spectral 
density  of  noise  iVo/2,  is  15  dB.  The  simulation  results  are 
drawn  from  4000  Monte  Carlo  runs.  The  symbols  of  each 
user  are  spread  by  Gold  codes  of  length  N  =  31.  The  first 
two  sequences,  mi  and  m2,  of  the  Gold  code  family  are 
the  primitive  m-sequences  with  octal  representation  of  their 
generator  polynomials  being  45,  75  and  their  loading  is  06, 
35,  respectively.  The  third  to  the  twenty-ninth  sequences  are 
generated  by  m»  =  mi  4-  Dx~z m2,  i  =  3,  33,  where 

+  denotes  bitwise  exclusive-or  operation  and  Dj  denotes  a 
phase  shift  of  the  m-sequence  by  j  units  [8],  The  weighting 
matrix  U  defined  in  (A.4)  was  employed.  Fig.  1(a)  and  (b) 
shows  the  bias,  -E[At]  (with  reverse  sign  to  be  shown  in 
log  scale),  and  the  mean-squared  error  (MSE),  E[Ar2],  with 
Gold  code  sequences  m2o  to  m24  and  delays  of  (10.67,  2.13, 
26.46,  27.00,  15.69,  16.21)  chips.  Each  interfering  user  is  10 
dB  stronger  than  the  desired  user.  The  scheme  in  Fig.  1(c) 
and  (d)  is  the  same  except  that  a  different  set  of  spreading 
codes,  mu  to  mi5,  and  delays  of  (25.40,  22.90,  26.95,  9.30, 
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(c)  (d) 

Fig.  1.  Approximation  and  simulation  of  estimator  performance  with  five  users:  (a)  and  (b)  are  for  the  first  set  of  code  sequences  and  delays,  and  (c)  and 
(d)  are  for  the  second  set  of  code  sequences  and  delays.  Dashed  line  with  o  is  simulation;  solid  line  is  approximation  by  (14)  or  (21);  dashed-dotted  line  and 
dotted  line  are  approximation  of  mean-squared  error  and  asymptotic  Cram&r-Rao  bound  by  Parkvall  [7],  respectively. 


8.28)  chips  are  used.  The  accuracy  of  the  approximation,  which 
depends  on  the  convergence  rate  of  the  Taylor  series  expansion 
to  the  true  J(r,  E5),  is  dependent  on  the  particular  realization 
of  spreading  codes  and  delays.  This  convergence  has  not  been 
fully  studied  at  this  time. 

Fig.  2(a)  and  (b)  shows  the  same  scheme  as  in  Fig.  1(a)  and 
(b)  with  five  additional  interfering  users,  each  with  energy  10 
dB  stronger  than  the  interested  user.  They  are  spread  by  the 
sequences  m25  to  m29  with  delays  of  (16.21,  19.67,  22.48, 
29.05,  11.31)  chips.  Discrepancy  between  the  approximations 
by  (14),  (21),  and  the  simulated  results  appears  to  be  far 
less  than  an  order  of  magnitude.  Fig.  2(c)  and  (d)  shows 
the  same  scheme  as  in  Fig.  2(a)  and  (b)  with  the  energy  of 
each  interfering  user  20  dB  stronger  than  the  interested  user. 
The  approximation  by  (14),  while  incidentally  achieving  better 
accuracy  in  the  ten-user  channel  than  in  the  five-user  channel, 
is  perceivably  better  than  that  of  [7]  in  all  four  situations.  It  is 
observed  that  both  the  estimator  performance  and  its  approxi¬ 
mation  are  fairly  insensitive  to  the  energy  variation.  However, 
they  are  sensitive  to  the  change  of  code  sequences  and  delays. 

V.  Conclusions 

In  this  paper,  we  have  provided  further  analysis  of  the 
statistics  of  the  modified  MUSIC  estimator  [3];  in  particular, 
approximations  of  the  first  and  second  moments  of  the  esti¬ 
mation  error  are  developed.  The  study  was  motivated  by  the 
need  for  delay  estimator  statistics  to  drive  the  construction 
of  delay-error  insensitive  MMSE-based  multiuser  receivers 


[6].  The  results  have  been  compared  to  a  previously  derived 
approximation  [7]  and  simulation  data.  In  addition  to  the  bias 
approximation,  the  new  approach  appears  to  provide  improved 
accuracy  of  the  approximation  of  the  mean-squared  error. 

Appendix 

SJ 

A.  Derivation  of  the  Gradients  g*,  q*,  and  the  Matrices  H,  H 

For  simplicity,  we  drop  the  subscript  p  throughout  the 
derivation.  While  the  derivation  is  valid  for  all  users,  the 
indexes  used  in  the  derivation  are  for  user  1  to  simplify 
expressions.  The  matrix  A  is  defined  as  [ai  a2]  with  ax  and  a2 
given  by  (3)  and  (4)  when  k  =  1.  The  first  and  second  partial 
derivatives  of  A  with  respect  to  the  fractional  delay  6  are 

;  a  dA 
A  98 

=  ^[(D^-D^Cx  (d^-D^Jcx 

(AI) 

»  a  dA 
A=  96 

=  0.  (A2) 

With  the  subscripts  of  pi  and  cx  dropped,  the  squared  norm 
of  the  code  sequences  as  a  function  of  6  is  given  by 

a? ai  =  c* [|J(2I +  21)  +1  c  (A3) 
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Sample  size  (M)  Sample  size  (M) 


(c)  (d) 

Fig.  2.  Approximation  and  simulation  of  estimator  performance  with  ten  users:  the  interfering  users  are  10  dB  stronger  than  the  interested  user  in  (a) 
and  (b)  and  20  dB  in  (c)  and  (d).  Dashed  line  with  o  is  simulation;  solid  line  is  approximation  by  (14)  or  (21);  dashed-dotted  line  and  dotted  line  are 
approximation  of  mean-squared  error  and  asymptotic  Cram&r— Rao  bound  by  Par kv all  [7],  respectively. 


for  i  =  1  and  2,  where  I  is  the  identity  matrix 

D,  =  (D“)Vrl+K1M>)"Dj; 

and 

d2=  (d^)  VP!+1)  +  (d^+1))"d^. 

In  the  following  discussion,  we  adopt  the  weighting  matrix 
in  [3]: 


with  the  subscript  p  of  Up  dropped.  By  substituting  (A3)  into 
(A4),  followed  by  taking  derivatives  with  respect  to  S,  we 
obtain  the  following  quantities: 


d  1  "Tat**  -t u 

ter  m*ri?  <A3> 

U(M)  =  |U(M)=^^  (Afl 


U(t\  i)  =  i)  = 


frxjKafaQ  u[  -  uf] 
(a"  a  <)4 


for  i  =  1  and  2,  where 

•“=  (i  4)c"<2i-b')c 

<  =  ^(21-DOc. 


(A8) 

(A9) 


The  (2,  1)  and  (1,2)  elements  of  U  and  thus  their  derivatives 


are  zero. 

Next,  we  will  find  the  derivatives  of  the  cost  function 
Jp(tp,  E»)  with  respect  to  6.  With  subscript  p  dropped,  the 
cost  function  J  is  given  by  J(r,  e*)  =  tr{UAHE„E^A}  = 
tr{UAH  (I  —  ESE^)A}.  Employing  the  derivatives  of  A 
and  U  given  by  (Al),  (A2),  (A5),  and  (A6),  finding  the 
derivatives  of  J  is  straightforward 


j(r,  ea) 

a  dJ(rt  E.) 

~  as 

=  tr{fJAKEnE"A  +  U(AHenE*A  +  AHEnE"A)} 
=  0  (A10) 

J(r,  Eg) 

A  a  j(r,  Eg) 
as 

=  tr{uAHEnE"A  +  2u(AHEne"  A  +  AHEnE*  A) 
+  2UA"Ene*A} 

=  tr{2UAHEnE"A}  (All) 


f 


138 

J\t,  E.) 

a  OJ(t,  Es) 
d6 

=  tr{UA"E„E"A  +  3U(A//enE"  A  +  A*E„E"  a) 

+  6UAHenE"A} 

=  tr{6UA"EnE"A}.  (A  12) 

The  gradient  of  a  complex  scalar  function  /(v)  is  defined 
•  as  the  directional  derivative  along  the  direction  pointed  by  a 
vector,  say  v0.  Thus,  the  gradient  VVo/(v)  is  a  vector  whose 
zth  element  is  df/dv0(i).  As  part  of  the  system  of  complex 
vector  operations,  the  convention  related  to  the  derivatives  of 
a  complex  number  z,  dz/dz  =  1,  dzjdzH  —  dzH /dz  —  0 
[12],  is  applied.  By  taking  the  gradient  of  j(r.  Es)  and 
J(r,  Es)  along  the  signal-space  eigenvector  e*,  we  obtain  the 
expressions  for  g i  and  qt- 

g,  =  VeJ(r,  Es) 

=  ^-tr{lJAHEnE*A  +  UAHEnE"A 
+UA"EnE"A| 

=  -  -^trjlJA^EsEf  A  +  UA"E3Ef  A 

oEii  t 

+UAHE3EfA} 

=  -  {a*uAt  +  A*uAr}e; 

=  ¥ej  (A  13) 

q,  =  Ve,j‘(r,E,) 

=  -  ^-tr{uAHE,Ef  A 

+  2u(AHEsEf  A  +  AffEgEf  A) 

+  2UAHE,Ef  A} 

=  -  {a*uat  +  2^A*UAt  +  a*uAt 

+A*UAT)}e*.  (A  14) 

Next,  we  consider  the  derivatives  of  the  Hessian  matri¬ 
ces  in  (17).  The  (&,  Z)th  element  of  H ij  is  defined  by 
H ij(k,  l )  =  d/de*(k)(dJ/ej(l))  [10],  where  the  inner  par¬ 
tial  derivative  is  equal  to  the  Zth  element  of  (A  13)  for  .  Since 
dJ/dej(l)  is  only  a  function  of  e^,  d/de*(k)(dJ/dej(l))  is 
\£(Z,  h)  for  t  =  j  and  zero  otherwise.  For  the  same  reason, 

l)  =  d/deiWid/djejil))  =  0,  VJk,  l  because  the 
expression  inside  the  parenthesis  is  a  function  of  only. 
Therefore,  H  and  H  are  given  by  (19)  and  (20). 
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Abstract -  The  problem  of  multi-user  detection  in 
Direct- Sequence  Code-Division  Multiple- Access  sys¬ 
tems  with  channel  state  mismatch  is  addressed.  An 
approximate  maximum-likelihood  criterion  is  derived 
based  on  knowledge  of  the  statistical  properties  of  the 
residual  channel  estimation  errors.  The  time  of  arrival 
is  assumed  to  be  a  constant  during  the  time  frame  of 
interest.  Errors  in  both  timing  information  and  fading 
coefficient  information  are  assumed.  A  Viterbi-like  re¬ 
cursive  demodulation  scheme  is  derived  and  analyzed. 
Performance  of  the  new  algorithm  is  studied  via  sim¬ 
ulation  and  the  development  of  performance  bounds. 

I.  INTRODUCTION 

Multi-user  detection  for  uncertain  environments  has 
evolved  in  two  directions.  The  first  approach  is  based 
on  joint  data  detection  and  parameter  estimation,  which 
is  computationally  challenging  especially  when  nonlinear 
parameters,  such  as  time  delays  are  to  be  estimated  [1], 
When  bandwidth  efficient  chip  pulses  are  employed,  the 
sampled  and  filtered  observation  signals  are  non-linear 
functions  of  the  relative  multi-user  time  delays.  Thus  we 
term  such  parameters  as  nonlinear  parameters.  Another 
approach  is  to  perform  blind  or  decision-directed  param¬ 
eter  estimation  followed  by  detection  algorithms  which 
are  optimized  under  the  assumption  of  perfect  parame¬ 
ter  estimates.  In  practice,  especially  in  dynamic  channels 
which  are  typically  experienced  in  mobile  communication 
systems,  perfect  estimation  can  hardly  be  achieved.  Mis¬ 
match  in  parameter  estimates  might  greatly  impact  detec¬ 
tor  performance  [2,  3].  In  this  paper,  we  pursue  modifica¬ 
tion  of  the  second  approach  to  reduce  detector  sensitivity 
to  parameter  mismatch. 

Due  to  the  nonlinear  nature  of  timing  in  detector  for¬ 
mulation,  isolation  of  the  timing  estimation  is  desired 
to  reduce  the  computational  effort.  Without  any  pre¬ 
sumptions  based  on  the  tracking  mechanism,  we  model 
the  residual  errors  in  the  timing  estimates  as  indepen¬ 
dent  Gaussian  random  variables.  The  flat  fading  process 
is  modeled  as  a  second-order  autoregressive  (AR)  pro¬ 
cess.  Efforts  towards  reducing  the  impact  of  imperfect 
parameter  estimation  have  been  previously  considered  in 
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Figure  1:  Receiver  diagram 


[3,  4,  5,  6].  Most  of  the  previous  works  for  asynchronous 
channels  considered  single  or  multi-shot  type  receivers 
where  edge  effects,  defined  as  a  partial  observation  of  the 
data  at  the  edge  of  the  observation  window,  contribute  to 
the  degradation  of  detection  performance.  For  additive 
white  Gaussian  noise  (AWGN)  channels,  we  proposed  a 
trellis-based  sequence  detector  based  on  an  approximate 
maximum-likelihood  (AML)  metric,  which  incorporates 
statistics  of  residual  timing  errors  [7].  In  this  paper,  we 
consider  extension  of  these  ideas  to  flat-fading  channels. 

This  paper  is  organized  as  follows.  The  system  model 
is  discussed  in  Section  II.  In  Section  III,  the  receiver 
structure  is  presented.  Bounds  on  performance  are  devel¬ 
oped  in  Section  IV  and  numerical  results  are  provided  in 
Section  V.  Final  conclusions  are  drawn  in  Section  VI. 

II.  SYSTEM  MODEL 

Processing  of  the  continuous-time  baseband  signal  is 
illustrated  in  Figure  1.  We  shall  assume  coherent  recep¬ 
tion  of  the  signal,  which  is  transmitted  through  channels 
without  resolvable  multipath  The  received  signal  is 
given  by 

K  M 

=  YL  ak(pi)bk(m)sk{t-mT-Tk)+n(t),  (1) 

Jk=l  m=—M 

where  a*(m)  is  the  amplitude,  6*  (m)  is  the  data  with  dif¬ 
ferential  BPSK  modulation,  r*  is  the  delay,  and  n(t)  is  an 

1  While  multipath  channel  models  are  more  appropriate  for  wide¬ 
band  signaling  such  as  DS-CDMA,  we  shall  study  flat  fading 
channels  as  a  necessary  precursor  to  the  study  of  multipath 
channels. 


AWGN  random  process.  Both  n(t)  and  a*(m)  are  com¬ 
plex  Gaussian  with  values  varying  in  time.  The  spreading 
waveform  sk(t)  is  given  by 

N—l 

sk(t)  =  ^2Mm-jTc),  (2) 

j=o 

where  a  chip-pulse  waveform,  ip(t),  is  modulated  by  the 
binary  spreading  sequence  0k(j)  The  square-root  raised 
cosine  pulse  shape  is  employed  for  ip(t).  The  front-end 
processor  (FEP)  at  the  receiving  end  is  a  linear  filter  with 
impulse  response  ip*  (~t).  For  the  convenience  of  analysis, 
we  will  approximate  the  fading  process  by  a  seconder der 
AR  process  [8], 

ak(t)  —  Pkia<k(t  -  1)  +  Pk2dk{t  -  2)  +  w(t),  (3) 

where  a*(t)  is  the  kth  user’s  fading  coefficient  at  time  t, 
p’s  are  the  parameters  in  the  AR  process,  and  w(t)  is  the 
Gaussian  process.  The  coefficients  in  (3)  are  given  by 

Pi  =  Ra(0)Jo(27T(i  -  l)t >T/cfc),  (4) 

where  Ra( 0)  is  the  average  power  of  the  process,  Jo(») 
is  the  zero-order  Bessel  function,  v  is  the  vehicle  speed, 
T  is  the  symbol  period,  c  is  the  speed  of  light,  and  fc 
is  the  center  frequency.  We  assume  that  the  fading  is 
slow  enough  so  that  ak  (m)  remains  a  constant  within  one 
symbol  period. 

The  vectors  of  delay  and  amplitude  are  written  as  the 
sum  of  their  estimates  and  residual  errors,  i. e. , 

r  =  f  4-  Sr 

a  =  a-f-5a,  (5) 


Figure  2:  Sampling  instants  at  the  output  of  FEP  is 
marked  by  the  asterisk.  ( qov  =  10;  qch  =  7). 

FEP  output  for  each  user  [7].  The  qsym  samples  for  each 
symbol  is  Tc/qov  apart,  where  qov  is  denoted  the  over- 
sampling  factor  with  qov  >  qsym •  The  sampled  signal  for 
a  channel  with  perfect  timing  and  amplitude  information 
may  be  expressed  as 

y  =  RAb  -I-  n,  (6) 

where  y  is  the  vector  of  CMF  outputs;  R  is  the  cross¬ 
correlation  matrix,  A  is  the  diagonal  amplitude  matrix,  n 
is  noise  vector  with  zero  mean  and  covariance  matrix  En, 
denoted  by  0,  En);  and  b  is  the  data  vector  (e.p.[10]). 
The  received  signal  is  arranged  in  the  order, 

y  =  [yf(-M) . . .  y £(-m)  . . .  yf (M) . . . y (7) 

III.  RECEIVER  STRUCTURE 

For  imperfect  timing  and  amplitude  information,  a 
linear  approximation  of  the  signal  portion  of  y  in  (6)  is 
obtained  via  a  Taylor  expansion, 


where  ST  and  5a  are  zero-mean  Gaussian  vectors  with  in¬ 
dependent  elements,  i.e.,  ST  ~  Af(0,Er),  <5a  ~  A/X0,£a) 
with  diagonal  covariance  matrices  ST  and  Ea.  This  scheme 
corresponds  to  adopting  the  maximum  average  uncer¬ 
tainty  principle  subject  to  fixed  variance  for  the  residual 
errors  [9]. 

For  bandwidth  limited  signals,  Nyquist  sampling  pro¬ 
vides  sufficient  statistics  for  data  detection  regardless  of 
the  time  of  arrival.  By  the  same  token,  Nyquist  sam¬ 
pling  takes  no  advantage  of  timing  information  when  it  is 
available.  If  the  timing  is  perfectly  known,  symbol-rate 
sampling  is  sufficient  for  data  detection.  In  this  paper,  we 
assume  knowledge  of  the  timing  estimates  prior  to  the  de¬ 
tection  procedure.  It  is  further  assumed  that  these  timing 
estimates  contain  residual  errors.  As  the  residual  errors 
are  small,  it  is  beneficial  to  slightly  modify  the  scheme  of 
symbol-rate  sampling.  The  proposed  sampling  scheme  is 
shown  in  Figure  2.  The  output  of  each  code  matched  filter 
(CMF)  is  sampled  non-uniformly  at  a  rate  of  qsym  times 
per  symbol,  or  equivalently,  q8ym  times  per  chip  at  the 


y  «  RBa  4*  (RB)<5a  -I-  (RBA)5r  +  n 

=  yo  +  Q0(b)(Ja  +  Qr(b)<Jr  +  n,  (8) 


where  B  =  diag(b).  The  matrix  R  is  obtained  by  substi¬ 
tuting  the  true  signal  space  after  sampling  with  the  signal 
space  corresponding  to  estimated  timing;  similarly,  sub¬ 
stituting  the  true  signal  space  with  its  derivative  with 
respect  to  timing  would  give  R.  That  is, 


R(r) 

R(r) 


®-(t)It=t 

0R(r)l 


dr 


I T-T 


(9) 


In  this  paper,  R  and  R  are  rectangular  matrices  since 
Qsym  >  1  is  Used. 

By  applying p(5a)  ~  ^(0,  Ea),  p(ST)  ~  jV(0,  ET)  and 
(8),  we  may  write  the  approximate  likelihood  function  as 


Pb(y)  ~  A/'Cyo,  E6), 


(10) 


<•* 


where  yo  =  RAb  and 

s6  =  Sn  +  Q0(b)E.Q?(b)  +  Qr(b)ErQf  (b).  (11) 
Thus  the  ML  estimate  of  b  is  given  by 
b  =  axg  maxfid(b) 

b 

=  axg  min  {(y  -  yo)*^  -  yo)  +  log  |S6|} . 

(12) 

The  cost  function  for  the  original  ML  formulation  with 
perfect  knowledge  will  be  denoted  by  ftc(b).  As  we  are 
considering  asynchronous  transmission,  the  complexity  of 
the  ML  detection  in  (12)  is  exponential  in  the  number  of 
total  symbols  transmitted.  It  is  noted  that  the  Viterbi  al¬ 
gorithm  (VA)  can  not  be  applied' directly  to  this  problem 
since  E&  in  (12)  is  sequence  dependent.  We  may  trans¬ 
form  the  received  signal  y  into  a  new  random  vector  z 
according  to  an  innovations-based  transformation  (IBT) 

in],  ,  s 

Z  =  r(y  -  yo},  (13) 

where  yo  is  the  estimated  mean  of  y  and  T  is  a  causal, 
whitening  transformation  matrix.  The  pdf  of  z  is  given 
by 

pb( z)  =  7r“(2M+1)^vm  exp{— zHz}.  (14) 

Since  the  elements  in  z  are  independent  and,  further¬ 
more,  they  are  not  explicitly  sequence  dependent,  a  con¬ 
ventional  VA  can  be  applied  in  terms  of  z  with  slight 
modification.  It  is  noted  in  [10]  that  z  is  causal  function 
of  n  and  y ,  but  not  of  b.  However,  we  may  still  apply  the 
VA  with  future  K  symbols,  as  well  as  the  past  K  symbols, 
represented  in  the  states  at  each  node.  The  trellis  for  a  2- 
user  system  is  shown  in  Figure  3.  The  size  of  the  trellis  is 
equal  to  that  for  a  channel  with  inter-symbol  interference 
memory  of  2 K  symbols. 

IV.  PERFORMANCE  ANALYSIS 

In  the  following,  the  detector  based  on  the  original  VA 
is  denoted  by  conventional  maximum-likelihood  sequence 
detector  (C-MLSD);  and  the  modified  version  presented 
in  Section  III  is  by  modified  MLSD  (M-MLSD).  Pairwise 
error  probabilities  (PEP’s)  are  obtained  by  considering 
the  situation  when  the  detector  erroneously  chooses  b 
instead  of  the  transmitted  data  b.  For  a  fixed  realization 
of  delay  errors,  we  will  present  the  PEP  only  for  the  C- 
MLSD  in  Rayleigh  fading  channels.  Comparable  analysis 
for  the  C-MLSD  in  AWGN  channels  can  be  found  in  [2]. 
The  difficulty  in  calculating  the  PEP  for  the  M-MLSD 
stems  from  the  following:  the  PEP  is  given  by  the  average 
probability  of  the  following  event  over  the  pdf  of  a  and  Saj 

nd(b)  -  Oi(b)  >  0,  (15) 


[bf<0),b,(1)*,(1).b,<2)]  [b1(1),b,(2>.bI(2),b1<3)] 

lb,(1).b,(1).b  ,<2)^(2)] 


Figure  3:  Top:  transmitted  data  sequence  by  two  users; 
Bottom:  corresponding  trellis  diagram  of  the  MLSD. 

where  fld(b)  is  given  in  (12).  Since  E&  is  a  quadratic  form 
in  a,  where  a  is  a  Gaussian  vector,  the  overall  expression 
in  (15)  is  a  complicated  function  of  a,  which  can  not  be 
analyzed  by  conventional  second-order  techniques. 

By  substituting  y  =  RB(a  +  <Sa)  +  n  and  yo  =  RBa, 
we  may  write  the  pairwise  error  event  as 

0  <  Ccr  =  ftc(b)-nc(b) 

=  X^PfcX  -  X^PjX  =  X^PcrX  (16) 


where 


X  = 

>t  si  »r]r 

\H  , 

s-1 

(R-R)B  RB  i] 

P  6  = 

(R  -  R)B  RB  I 

.  m 

H  1 1 
s-1 

»  _  j 

Pi  = 

RB  -  RB  RB  I 

RB-RB  RB  ij 

(17) 

All  the  quantities  in  (17)  jure  real,  except  x  which  is  a 


(18) 


complex  Gaussian  vector.  The  PEP  is  given  by 

*%.>!!]=  e  n  ('-¥) 

i  3  '  ' 

>  o  i  » 

where  Aj’s  axe  the  eigenvalues  of  £xPcr,  with  £x  being 
the  covariance  matrix  of  x  [6].  It  is  assumed  in  (18)  that 
none  of  the  eigenvalues  are  identical. 

A  practical  example  for  Sx  can  be  obtained  by  con¬ 
sidering  approximation  of  the  flat  Rayleigh  fading  pro¬ 
cess  in  (3).  We  assume  an  ideal  scenario  where  a*(ro)  = 
pk\dk{m-  1)  -bp*2a*(m-2)  and  5a,k{m)  =  e(t);  that  is, 
the  estimator  is  capable  of  tracking  the  predictable  part 
of  the  fading  coefficients.  Thus,  the  corresponding  PEP 
is  a  lower  bound  of  that  for  a  practical  system  which  em¬ 
ploys  a  causal  channel  tracking  process.  According  to  the 
AR  model,  6a,k (m)  is  a  white  Gaussian  process  with  co- 
variance  matrix  diag([crf , . . . ,  <7^,  af ,  •  *  * ,  cr^, . . .]),  where 
o\  —  £{|Jat*(m)|2}.  We  assume  that  the  fading  process  is 
independent  amongst  the  users.  Thus,  we  focus  on  one  of 
the  users  and  drop  the  index  k.  The  correlation  function 
of  d(m)  =  a(m)  —  Sa(m)  is  given  by 


Figure  4:  BER  and  analytical  bounds  vs.  SNR  for  C- 
MLSD  and  M-MLSD  in  flat  Rayleigh  fading  channels. 

These  bounds  are  accurate  for  a  sequence  detector  only 
when  data  of  infinite  horizon  is  considered.  To  reduce 
the  computational  effort,  a  finite  horizon  of  five  data  bits 
(M  =  2)  per  user  is  considered  with  the  interested  bit, 
&i(0),  in  the  middle.  Thus,  the  bounds  presented  here 
are  expected  to  be  slightly  higher  than  their  true  values, 
which  correspond  to  the  infinite  horizon  case. 


Ra(n)  =  E{a(n)a*(n  -  m)} 

=  E  {[pia(m  -  1)  +  p2a(m  -  2)] 

\pia{m  -  1)  +  p2a(m  -  2)]*} 

=  (IPil2  +  \p2\*)Ra(rn)  +p1p*2Ra(m  +  1) 
+p\p2Ra{m- 1).  (19) 


Given  that  E„  and  Sj0  are  known  and  E„  can  be  ob¬ 
tained  from  (19),  the  covariance  matrix  of  x  is  given  by 


0 

0 

Sn 


(20) 


To  determine  a  lower  bound  of  performance,  consider 
a  scheme  where  the  vector  b  is  transmitted.  Assume  a 
genie  tells  the  receiver  that  either  b  or  b m  is  transmitted, 
where  i>\{  =  arg  max  P{fid(b)  —  fi<*(b)  >  0|b}.  That 

b 

is,  bM  is  the  most  probable  symbol  to  be  mistakenly  de¬ 
tected  when  b  is  transmitted.  The  BER  is  bounded  by 
the  performance  of  this  genie-aided  receiver, 

BER>  >  2(^+1)*  ^  "t*  Ppw' 

b  &1(0)S1(i)#  MOjbjd) 


where 

ppw  =  p[nd(b)  -  nd(b)  >  o|b].  (21) 

The  upper  bound  employed  is  simply  the  union  bound, 

BERl  <  2 (2A?TW  Ppw' 

b  b 

h(O)hd) 


V.  NUMERICAL  RESULTS 

In  this  section,  preliminary  results  are  obtained  under 
the  assumption  that  the  only  source  of  prediction  .error 
in  the  fading  coefficients  is  the  innovations  process  w(t) 
and  the  ambient  noise.  Thus  the  channel  estimator  per¬ 
fectly  tracks  the  predictable  portion  of  the  AR  process  in 
(3).  Therefore,  the  results  presented  here  provide  a  lower 
bound  of  the  bit-error  rate  (BER).  Due  to  complexity 
of  the  algorithm,  a  two-user  system  with  Gold  codes  of 
length  31  is  considered.  The  decision  depth  of  the  simu¬ 
lated  MLSDJs  is  10  times  the  channel  memory.  The  SNR 
is  defined  as  the  peak  signal  energy  at  the  CMF  output 
over  the  variance  of  the  noise.  Due  to  timing  uncertainty, 
the  “effective”  SNR  is  lower  than  that  with  regular  SNR 
definition. 

Figure  4  shows  the  performance  in  flat  Rayleigh  fad¬ 
ing  channels  approximated  by  a  second-order  AR  process 
with  the  following  system  parameters:  center  frequency 
fc  =  900  MHz,  vehicle  speed  v  =  100  km/hr,  symbol  rate 
T~l  =  9.6Ksec~l,  and  NFR=4  dB.  We  acknowledge  the 
fact  that  the  accuracy  of  timing  estimation  should  im¬ 
prove  with  increased  SNR,  and  thus  o\k  should  change 
with  SNR.  However,  without  studying  any  specific  tim¬ 
ing  estimators  in  this  paper,  we  consider  fixing  the  value 
of  o\  as  an  artificial  but  acceptable  choice.  The  sys¬ 
tem  is  tested  with  fixed  tracking  error  of  0.1TC,  which  is 
equal  to  the  standard  deviation  of  tracking  errors  used 
to  construct  the  receiver.  To  reduce  the  error  of  the  lin¬ 
ear  approximation  for  large  SNR,  we  found  the  param¬ 
eters  qaym  =  7  and  qov  =  10  to  be  sufficiently  large  for 


